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ABSTRACT 

A  prototype  coupled  Marina  Atmospheric  Boundary  Layer 
(MABL)  and  Oceanic  Boundary  Layer  (OBL)  model,  from  physical 
models  which  have  been  separately  formulated,  is  developed. 
The  observational  efforts  will  be  directed  toward  coupling 
(local)  point  MABL  and  OBL  features.  Emphasis  in  the  char- 
acterization of  local  features  will  be  on  the  evolutions  of 
the  adjacent  well-mixed  layers.  The  approach  will  be  to 
compare  observed  evolutions  in  the  oceanic  and  atmospheric 
boundary  layers  with  predictions  from  bulk  models  wherein 
the  ocean  influences  the  atmosphere  through  the  surface  tem- 
perature which,  in  turn,  is  influenced  by  atmospheric  forc- 
ing. The  atmospheric  forcing  parameters  are  the  local 
surface  ^urbulent  fluxes  and  radiation.  Surface  turbulent 
fluxes  are  controlled  by  the  surface  temperature  and  atmos- 
pheric mixed  layer  parameters.  Radiation  is  con-rolled  by 
predicted    cloud    formation    and    dissipation. 
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I.        INTRODUCTION 

The  idea  that  both  atmospheric  and  subsurface  environ- 
mental conditions  can  have  a  significant  effect  (enhancement 
or  degradation)  on  modern  weapon  and  sensor  systems  has  now 
become  a  reality  to  many  military  leaders.  This  has  led  to 
a  greater  demand  for  more  reliable  and  accurate  atmospheric 
and  oceanic  prediction  schemes.  Prediction  of  the  existence 
and  evolution  of  electromagnetic  ducts  are  important  in  pre- 
dicting over  the  horizon  surface  to  surface  detection.  Pre- 
dictions of  the  inversion  height  is  important  where  it  is 
known  that  optical  propagation  is  degraded  due  to  turbulence 
and  extinction.  Ocean  mixed  layer  predictions  are  important 
for  determining  convergence  zone  distances,  optimum  source 
level  depth,  and  cutoff  frequency  for  acoustic  surface  duct 
propagation. 

The  ability  to  provide  accurate  forecasts  requires  a 
complete  understanding  of  the  dynamical  processes  occurring 
in  each  of  the  atmospheric  and  oceanic  boundary  layers,  and 
as  important,  the  processes  associated  with  coupled  local 
changes    in    the   two   layers. 

Very  little  is  known  about  the  interactive  nature  of  the 
adjacent    turbulent   boundary  layers.  The    objective    of   this 
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thesis  is  to  formulate  a  model  which  couples  two  already 
existing  and  proven  models  for  the  respective  adjacent 
well-mixed  boundary  layers.  This  coupling  is  achieved  by 
matching  the  surface  boundary  conditions  for  stress,  mass 
flux   and   hea-1-    flux    (Figure   1)  . 

This  coupled  model  may  then  be  used  to  examine  the  pro- 
cesses of  coupled  local  oceanic  and  atmosoheric  mixed  layers 
at  time  and  spatial  scales  of  12-18  hoars  and  50-100  km, 
respectively.  Significant  changes  in  both  mixed  layers  can 
occur  at  these  scales.  Explaining  these  changes  requires 
coupled  evaluations  of  exis-ing  atmospheric  and  oceanic 
well-mixed    layers   at    a  point. 

At  the  present  time,  we  only  have  qualitative  evidence 
-that    changes      in   each   layer     are    correlated.  The   rigorous 

specifications  or  descriptions  of  the  coupling  will  require 
considerable  interpretive  efforts  and,  perhaps,  improved 
formulation   of   the   separa-e   models. 
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II.       BACKGROUND 

A.        GENERAL    BOUNDARY    LAYER    FEATURES 

The  Marine  Atmospheric  Boundary  Layer  (MABL)  extends 
from  the  surface  through  the  capping  inversion.  The  Oceanic 
Boundary  Layer  (OBL)  extends  from  the  sirface  to  the  sea- 
sonal thermociine.  The  air-ssa  interface  is  bordered  by 
oceanic  and  atmospheric  turbulent  mixed  layers  which  effec- 
tively insulate  the  quasi-geos trophic  oceanic  and  atmos- 
pheric flows  below  the  thermociine  and  above  the  inversion, 
respectively.  The  primary  source  of  the  turbulence  is  the 
velocity  (current)  and  buoyancy  (density)  gradients  created 
by  interaction  a+  the  interface.  The  large  vertical  mixing 
yields  nearly  homogeneous  (well- mixed)  wind,  "-emperat  ur  e  and 
humidity  profiles  above  the  surface.  In  the  ocean,  the  term 
mixed  layer  refers  to  the  region  cf  vertical  uniformity  in 
mean  velocity  and  density.  Although  density  is  a  function 
of  both  -emperature  and  salinity,  it  is  usually  dependent 
upon  temperature  except  when  large  near-surface  salinity 
changes  occur  due  to  heavy  precipitation  :r  ice  melt  in  high 
latitudes.         At    the      bottom   of    the   ocean    mixed      layer    and    at 
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the  top  of  the  atmosphere  mixed  layer  are  thin  transition 
regions  (thermocline  or  inversion)  .  The  mixed  layers  inter- 
act with  the  upper  and  lower  layers  at  these  regions  by 
means  of  turbulent  forced  entrainment.  Typical  model  pro- 
files depicting  the  boundary  layer  regions  for  both  the 
atmosphere   and  ocean    are    shown   in    Figure   2. 

The  well-mixed  nature  of  both  the  turbulent  MABL  and  OBL 
implies  features  for  the  vertical  distribution  of  mean  val- 
ues and  vertical  fluxes  of  wind,  temperature  and  humidity  in 
the  atmosphere  and  density  and  momentum  in  the  ocean.  One 
feature  is  that  the  properties  which  are  conserved  during 
mixing  can  be  treated  as  being  constant  with  height  (depth) 
in  the  boundary  layer.  A  second  feature  is  that  vertical 
fluxes  of  the  well-mixed  parameters  vary  linearly  with 
height.  These   features      enable      predictions      to    be      based 

solely  on  the  surface  and  inversion  (entrainment)  fluxes  if 
advection    is    negligible. 

Another  important  source  of  energy,  in  addition  to  tur- 
bulent kinetic  energy,  in  both  mixed  layers  is  radiation. 
Unlike  the  MABL,  most  of  the  solar  radiation  does  not  pene- 
trate the  OBL.  Therefore,  downward  turbulent  heat  flux  in 
the    OBL   is   as    important    as   the    upward    flux.         Since   the    long 
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Figure  2.   Simplified  Atmospheric  aid  Oceanic  Boundary  Lay< 
Temperature  Profiles. 
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wave  (infrared)  radiation  flax  is  absorbed  and  emitted 
within  the  first  millimeter  of  the  sea  surface,  this  flux 
component  is  essentially  another  component  of  the  oceanic 
surface  heat  flux.  The  total  surface  heat  flux  is  comprised 
of  sensible  and  latent  heat  and  the  long  and  short  wave 
solar  radiation.  Beth  the  long  and  short  wave  (solar) 
fluxes  into  the  ocean  are  critically  dependent  upon  the 
existance   of   clouds   in  the  atmosphere. 

In  both  the  HABL  and  OBL,  buoyantly-created  velocity 
fluctuations  and  mixing  have  the  most:  pronounced  effect  on 
mixed  layer  characteristics.  Garwood  (1977)  states  that 
buoyancy-forced  vertical  mixing  has  a  more  obvious  and 
direct  role  in  the  mechanical  energy  budge-,  as  shown  in 
Fiqure  3.  In  the  atmosphere,  even  under  near  neutral  condi- 
tions, buoyancy  forced  fluctuations  and  nixing  can  be  quite 
larqe.  These    buoyancy      driven    energetio      eddies    fill      the 

atmospheric  boundary  layer  from  the  surface  to  the  inversion 
where  they  entrain  warm,  dry  air  and  bring  more  momentum 
into  the  mixed  layer.  If  eitrainment  causes  the  MABL  to 
extend  above  the  lifting  condensation  level,  clouds  or  fog 
will   form    within  the    layer. 
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Figure  3.   Mechanical  Energy  3udget  for  the  Ocean  Mixed 
Layer. 

It  is   the  realization  that  buoyancy   driven  entrainment 

depends  on   the  near   surface  stability,    which  in   turn  is 

controlled  by  local   atmospheric  forcing  Df  the   upper  ocean 

layer  which   makes  a  coupled  pea  diet  ion   approach  necessary. 

A  circularity   in  cause   and  effect   arises  because   surface 

buoyancy   flux  induced  entrainment  not   only  increases   the 

depth  of   the  JIABL,   but  it   also  changes  its  effect   on  the 

ocean  mixed  layer.     The  role  of  the   surface  buoyancy  flux 

means  changes  in   the   oceanic   well-mixed  layer  could   be 

either  the   cause  or   the  result   of  entrainment.     Another 


18 


cause  and  effect  relationship  is  the  forma-ion  and  dissipa- 
tion of  clouds  within  the  ABL.  Clouds  can  be  caused  by 
changes  of  the  ocean  surface  temperature.  In  turn,  clouds 
have  a  profound  effect  on  the  short  and  long  wave  radiation 
budget   for   the   OBL. 

These  cause  and  effect  relationships  shew  that  the 
assessment  of  relative  roles  of  dynamic  processes  in  both 
layers  reguires  coupled  evaluations  using  relatively  complex 
physical  models.  A  data  set  collected  (Davidson  et  al  1983) 
over  a  48  hour  period  in  May  197  8  demonstrates  the  nature  of 
coupled  local  changes  in  the  atmospheric  and  oceanic  bound- 
ary layers.  For  a  period  from  1700  to  0500  local  time, 
clouds  existed  in  the  MABL  and  winds  increased  steadily  from 
6  to  9  ms~  during  this  period.  It  is  suggested  this 
increase  was  associated  with  deepening  of  the  MABL  (and 
therefore,  the  entrainment  of  additional  momentum  into  the 
mixed  layer)  which  is,  in  turn,  associated  with  the  surface 
temperature.  The  oceanic  boundary  layer  deepened  as  the 
wind  reached  a  maximum,  resulting  in  the  surface  temperature 
dropping    from    14.5°    to   13.0°C. 

There    were      other   atmospheric   response    and      forcing    fea- 
tures    observed      during      the      entire      48      hour      period      from 
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5/19/1200    to    5/21/1200      PST  which   are    believed      tc   have    been 
associated    with     coupling    between     the   two      layers.         Atmos- 
pheric and  oceanic  mixed   layer  observations    appear   in    Figure 
4. 
These  additional   features   are: 

a)  The    ABL   depth    increased   from    250    to    750    meters    over 
the   48    hour   period.    The   changes   occured    in    relatively 
short   intervals,    from    20/0000    to    23/0400,    and   from 
21/0000   to    21/0500.      Ths    ABL   depth    remained  at    near 
constant    depths  during    the   intervening    12    to    18    hour 
periods. 

b)  The   surface  layer   temperature    increase    from   20/0000 
to    20/1500   is    indicative    of   entrainment   of  overlying 
warm   air.      This  en  trains  e  n-1-    presumably    was    a    factor 
in    the    wind   speed    increase   over   the    same    period.    The 
speed    increase,   in   turn,     presumably    caused  the    warm 
shallow    OEL   to    be    destroyed  which    dropped  the   surface 
temperature   after    20/1800. 

c)  The    most    dramatic    change    in    the    flABL    occured   during 
the    period    from  21/0000   to    21/0500    which    was    immedi- 
ately  after   the  surface   temperature    had    dropped    due 
to    wind    mixing. 
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Figure  4.       A'mosDheric  and  Oceanic   Mixed    Layar   Observations 
during   CEWCOM-78.       (a)    pD~ential    ^emDeraturs 
profixes    ( I)    acoustic    sounder    raccrd"  (c)     10    mete: 
wind    speed    (d)     10    mater   and   surface   -tempera-1: ure 
and    (e)    XBT  traces. 
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Questions  which  cculd  be  posed  on  the  basis  of  the  above 
data   are: 

a)  Did  the    dramatic    change   in   The    ABL    depth    at   21/0000 
occur   as   a   result    of  the    ocean    forcing  the  atmosphere? 

b)  Does    the    fact    that    the    mixed    layer    evolved  with    rapid 
changes    follow€d    by  equilibrium   periods    indicate    feed- 
back  between    the    boundaries? 

c)  What    wculd   have   occurred    if   the    wind    had    reached    its 
maximum   value    at    a   time   other    than    after    sunset? 

d)  Was   the    time   of  occurrence   of    the    wind   maximum    and, 
hence,    the  ocean    mixing   ccn-rolled    by    the   coupled 
mixed    layer   system? 

This  thesis  will  attempt  to  answer  some  of  these  ques- 
tions on  the  basis  of  the  formulated  coupled  model.  As 
stated  earlier,  all  answers  can  only  bs  obtained  through 
considerable   interpretive   efforts. 

The  scope  of  the  overall  scientific  effort  is  illus- 
trated in  Figure  5.  The  local  oceanic  descriptions  will  be 
made  by  a  bulk  oceanic  mixed  layer  model  (the  Garwood 
Model).  The  atmospheric  descriptions  will  be  made  by  the 
Naval  Postgraduate  School  (NPS)  Marine  Atmospheric  Boundary 
Layer   model    (Davidson   et      al    1983).         Tactically    significant 
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operational  descriptions  (forecasts)  to  which  the  improved 
understandings  would  apply,  appear  on  the  extreme  right  hand 
side   of  Figure   5. 

B.        MARINE    ATMOSPHERIC    BOUNDARY    LAYER    (MABL)     MODEL 

The  NPS  boundary  layer  modal  is  a  zero-order,  two  layer, 
integrated  mixed  layer  model.  The  two  layers  consist  of  the 
well-mixed,  turbulent  boundary  layer,  and  the  relatively 
non-turbulent      free      atmosphere.  The   iDdel      is      based      on 

entrainment  energetics  formulated  by  Stage  and  Businger 
(1981)  and  en  radiative  transfers  described  by  Davidson  et 
al     (1983). 

The  inversion  (or  transition  zone)  separates  one  layer 
from    The    other.  In    a    zero-order   model,         the    inversion    is 

assumed  to  have  zero  thickness,  hence,  a  discontinuity,  or 
jump  occurs  in  the  profiles  of  the  conservative  parameters 
at    the  inversion. 

The  present  model  requires  as  input,  an  initial  atmos- 
pheric sounding,  the  mean  winds  at  a  level  within  the  sur- 
face layer  (10-30  meters)  and  the  surface  temperature. 
Boundary  layer  wind  values  for  the  duration  of  the  forecast 
period  can  be  obtained  from  larger  scale  models.  In  this 
study,      actual    hourly    wind   speed      measurements    were   input    as 
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part  of  the  model  initialization.  Sea-surface  temperature 
remains  unchanged  from  the  initial  input  value.  However,  in 
the  coupled  model,  the  sea-surface  temperature  will  change 
during  the  forecast  period  and  will  be  that  predicted  by  the 
OBL  model.  Well-mixed  temperature  and  humidity  are  predic- 
ted so  the  surface  wird  and  wind  shear  at  the  inversion  are 
the  only  atmospheric  variables  which  have  to  be  prescribed. 
The  larger  scale  subsidence,  obtained  from  the  large  scale 
model  must   also    be   prescribed   for   the    forecast   period. 

Three  methods  can  be  used  to  compute  the  large  scale 
vertical  velocity  (subsidence)  from  single  station  observa- 
tions. These  three  methods  are  well  described  by  Gieason 
(1982)  and  include  the  kinematic  and  adiabatic  methods  and 
the  integration  of  the  moisture  budget  equation  (Q-methcd) . 
Gieason' s  study  shewed  the  Q-method  displayed  the  most  merit 
as  a  si.i  gle-stat  ion  assessment  of  subsidence.  This  method 
was  used  to  compute  subsidence  for  this  study.  The  solar 
zenith  angle  used  in  the  computation  of  the  short  wave  radi- 
ative flux  is  determined  from  the  latitude,  Julian  day  and 
time   of   day  of   the    initial  input. 

The  model  predicts  the  time  evolution  of  the  inversion 
height,    and   the    mixed   layer   values   of    9    (temperature)       and    q 
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(humidity),  as  well  as  their  jumps,  at  3D-minute  intervals. 
The  prediction  model  is  run  out  to  24  hours.  The  steps  in 
the  prediction  computation  are  shown  in  Figure  6,  where  it 
is  noted  tha^  procedures  are  the  same  for  clear  and  cloudy 
cases  except  for  entrainment  computation  and  estimation  of 
cloud   top   cooling. 

The  model  uses  the  standard  integrated  rate  aquations 
(Ter.nekes  and  Dreidonks,  1981)  to  predict  the  time  evolution 
of  the  conservative  variables  and  their  jumps  at  the  inver- 
sion.     These    equations  are: 

h(Dx/Dt)=    (w'x')q    -(w'x'l     +    source  (1) 

h(D4x/Dt)=    hrx  (dh/dt)     -    (w«x»)0+     (v'x')h-    source       (2) 

where  C  is  the  lapse  rate  above  the  inversion  and  the  source 
term  is  equal  to  -  (Fnh-  Fno)  /fcp  for  x  =  9  and  equal  to  zero  for 
X=q.  Fn  is  the  net  radiative  flux.  The  subscripts  zero  and 
h  refer  to  surface  and  inversion  height  values,  respec- 
tively. 

Stage  and  3usinger*s  (1981)  entrainment  velocity  parame- 
terization is  used  to  close  the  system  of  equations  and 
determine  the  time  evolution  of  the  inversion  height.  The 
closure    assumption   states      that    the   dissipation    rate      of   the 
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TIME 
STEP  30 
JilN. 


NO 


COMPUTE  : 
ENTRAPMENT 


INITIAL  RADIOSONDE 

TEMPERATURE,  HUMIDITY 

SURFACE  TO  3  KM 


PRESCRIBE: 

SURFACE  WINDS 

SEA  SURFACE  TEMPERATURE 

SUBSIDENCE  (DIVERGENCE) 

ADVECTION 


COMPUTE: 

SURFACE  FLUXES 
CONDENSATION  LEVELS 


COMPUTE  MEW: 

MIXED  LAYER  DEPTH 
WELL  MIXED  TEMPERATURE 

AND  HUMIDITY 
JUMP  STRENGTHS 


COMPUTE  NEW: 
EM  DUCTS 
OPTICAL  TURBULENCE 

AMD  EXTINCTION 
DISPERSION 


YES 


COMPUTE : 
ENTRAPMENT 
CLOUD  TOP  COOLING 


Figure  6.      Schematic    cf   Input,    Prescription    and   Computing 
Steps   in    HABL    Prediction. 
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turbulent   kinetic  energy    (TKE)       is      a    fraction     (1-A)       of   the 
preduction    rate.         ■ A •   is   the    entrainment    coefficient,    taken 

as    0.2. 

Mathematical  modeling  of  radiative  flax  transfer  is  an 
extremely  complicated  but  still  developed  subject.  Unfortu- 
nately, even  the  simple  radiative  transfer  models  are  still 
extremely  complex  in  comparison  to  most  other  physical 
parameterizations  used  in  the  mixed  layer  model.  Uncer- 
tainty in  background  aerosols,  atmospheric  absorbing  gases 
(water  vapor,  carbon  dioxide,  ozone)  and  cloud  droplet  size 
spectra  are  sources  of  error  in  radiation  calculations. 
Long  and  short  wave  radiation  flaxes  are  computed  sepa- 
rately. 

The  long  wave  radiation  flax  calculation  was  modified  to 
permit  non-black  stratus  clouds  by  introducing  the  cloud 
emissivity,  £  ,  which  is  a  function  of  the  total  cloud  liq- 
uid content,  W.  Since  the  cloud  liquid  water  profiles  are 
approximately  linear  with  height  (Davidson  et  al,  1982), 
cloud  liquid  water  content  and  cloud  emissivity  are  given 
by: 

W=    0.5j?a(h    -   &c)q(h  (3) 

£c=    1    -    exp(-aW)  (4) 
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where  pa  is  the  density  of  air  (  1.25X  10"3  gm  cm"3  )  ,  h  is  the 
height  of  the  mixed  layer  (cloud  top),  Zc  is  the  lifting 
condensation  level  (cloud  bottom),  a=0.153  (Slingo,  1982), 
and    g)h     is   the    liquid    water   content    at    ths    cloud   top. 

The  long  wave  cloud  top  net  radiation  flux,  Lnh ,  can  be 
calculated  from  xhe  cloud  top  temperature,  Th ,  and  the 
effective  radiative  sky  temperature,  Tsky ,  using  the  Stefan- 
Boltzmann  Law.  At  the  cloud  bottom,  this  flux,  Lnc,  is  cal- 
culated     in        the     same        manner      using        the      sea        surface 


.emperature,   Ts  ,      and  cloud   bottom   temperature, 


fluxes  are   given   by: 


-c 


These 


Lnh-£c<T(Th     "    Tsky  ) 
Lnc  =  £c<T(Ts4-    Tc4  ) 


(5) 
(6) 


where  CT  is    Stefan's   constant    (U.61    X      10"11   )       and  £0  ,      cloud 
emissivity,    is   obtained    from    equation    4.         The    net    long    wave 

,    becomes: 


flux    at    the   surface. 


ong 


F.ong    =    tf-(?s4  "    £CT4-     (1-£c)Tsky4) 


(7) 


where  T   is      the   average    temperature   of   the      cloud.         For   the 
cloud   free   case,        the  net    fluxes   are    calculated      at    Z=h   and 
Z=0    using   the    standard  method      (Fleagle    and    Businger,       1978) 
of    integrating      the    flux   emissivity     profile.      The     net    long 
wave    flux    at   the   surface    for   the    clear   sky    case    is: 
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Flong    "    Fu     "    Fd  (8) 

where  Fu  is  the  upward  radiativa  flux  and  Fd  is  the  downward 
radiative    flux. 

Surface  flux  parameterization  (bulk  aarodynamicai)  for- 
mulae 

1/2 
U*=     Cd    U10  <9a> 

T*=    C(,1/2(0O-    9)  (9b) 

q*=    Cfl|/2(q0-    q)  (9c) 

are  used  to  determine  the  surface  fluxes  3f  momentum,  sensi- 
ble   heat    and   lament    heat.      Thesa    fluxes   are    given    by: 

2 

u'w's   u*  (momentum)  (10a) 

T'w'=   u*T*  (sensible    ha  at)  (10b) 

g»w«  =   u*q*  (latent    heat)  (10c) 

where  Cd  and  Cg  are  the  ten  mater  stability  dependant  drag 
coefficient,  9  is  -hz  potential  tamperatura,  and  q  is  the 
specific  humidity.  The  subscript  zero  danotes  the  surface 
value. 

The  shor-  wave  radiativa  flux  is  calculated  using  the 
delta-Edingtcn  method  (Joseph,  1976)  .  This  was  added  to  the 
model  to      account    for      heating   of   the      mixad    layer      by    solar 
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radiation.  Shor-1,  wave  extinction  is  dominated  by  scattering 
(as  opposed  tc  absorption).  This  scattered  radiation  forms 
a  second  short  wave  radiative  component  usually  referred  to 
as  diffuse'  solar  radiation.  Scattering  in  the  mixed  layer 
is  due  to  atmospheric  particles,  cloud  droplets  in  the 
cloudy  case  and  sea-salt  aerosols  in  the  cloud-free  case. 
This  MABL  model  computes  both  the  direct,  and  diffuse  radia- 
tive components  to  determine  the  total  short  wave  radiation 
flux  at  the  surface.  An  explanation  of  the  delta-Edington 
method  and  all  of  the  parameters,  atmospheric  factors,  and 
equations  involved  is  a  very  complex  subject  and  is  beyond 
the  scope  of  this  paper.  An  excellent  review  of  this  sub- 
ject has  beer,  published  by  Fairall  (1981).  An  important 
short  wave  radia-ive  parameter  that  is  prescribed  in  the  A3L 
model  is  the  fraction  of  reflected  short  wave  radiation,  Ag, 
from   the    sea    surface.      It    is   taken   as    0.1. 

C.       OCEANIC    BOUNDARY     LAYER    (OBL>      MODEL 

Garwood  (1977)  developed  an  ocean  mixed  layer  model 
using  the  Navier-Stokes  equation  cf  motion  with  the  geostro- 
phic  component  eliminated,  the  continuity  equation  in  incom- 
pressible water,  the  heat  equation  fron  the  first  law  of 
thermodynamics,  the  conservation  of  salt  equation,  and  a 
linearized   equation   of  state. 
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Prediction  of  the  rate  of  deepening  (d:  retrsat)  cf  the 
mixed  layer  is  dependent  upon  an  understanding  of  the  dynam- 
ics of  the  entrainmen4-  process.  The  turbalence  of  the  cver- 
lyinq  mixed  layer  provides  the  energy  needed  to  destabilize 
and  erode  the  underlying  stable  water  mass  (Garwood  1977). 
Therefore,  -he  turbulent  kinetic  energy  budget  is  ths  basis 
for  the  entrainment.  This  system  is  closed  using  a  mean- 
turbulent-field  modeling  of  the  vertically  integrated  equa- 
tions for  the  individual  TKE  components,  plus  the  inclusion 
of    the   bulk   buoyancy    and    momentum   equations. 

Seperate  ver-ical  and  horizontal  equations  for  TK2  are 
used  to  better  handle  the  mixing  process.  Buoyancy  produced 
energy  is  somewhat  more  efficient  than  shaar  production  as  a 
source  of  energy  for  vertical  mixing  because  of  its  unique 
effect  on  the  vertical  component  of  the  turbulent  velocity. 
The  buoyancy  equation  is  generated  from  the  heat  and  salt 
equations    together   with   an  equation   of   state: 

}  =  _fo  [l-«(9    -    9„)    +   0  (S   -   Sa  )]  (11) 

and    buoyancy    is   given   by: 

b=    g(j>0   -/)//o  (12) 
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where  6  is  temperature,  S  is  salinity,  o  density,  g  gravity, 
and  the  cons-ants  <x  and  S  are  the  expansion  coefficients  for 
heat  and  salt,  respectively.  The  ■'•ilde  represents  the  total 
instantaneous  value  and  the  subscript  zero  denotes  an  arbi- 
trary, but  representative,  constant  value.  Using  b  instead 
of  9  (temperature)  for  the  definition  of  buoyancy  allows 
this  model  to  be  applicable  in  situations  where  evaporation 
and  precipitation  contribute  significantly  to  the  surface 
buoyancy  flux.  Temperature  is  in  most  oases  considered  to 
be  the  dominating  factor  in  the  density  profile.  The  rela- 
tive importance  of  the  salinity  effect  on  the  short  term 
evolution  of  the  density  profile  is  generally  nor  signifi- 
cant    (Miller,    1976). 

Temperature  and  salinity  profiles,  and  for  extended 
forecasts,  the  wind  driven  horizontal  current  profiles,  are 
reguired  as  input  for  model  initialization.  These  are  used 
to  compute  the  mixed  layer  depth.  The  Garwood  OBL  model 
defines  the  mixed  layer  depth,  h  ,  as  the  shallowest  dep'rh 
at  which  the  observed  densitv  value,  c7~{  ,  is  Q.02(Tt  units 
greater   than  the    observed    surface   densi-y    value. 

Ocean  environmental  parameters  that  have  ^o  be  pre- 
scribed     include;        the      fraction  of      short      wave      radiation 
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absorbed  in  the  top  one  meter  of  ■'-tie  ocean,  the  radiation 
extinction  coefficient  for  det erminatioi  of  the  downward 
turbulent  heat  flux,  and  the  critical  Richardson  number  for 
a   stability  adjustment  at   the   bottom    of   the    mixed   layer. 

The  parameters  required  for  surface  boundary  condition 
computations  include  wind  speed  and  direction,  cloud  cover, 
sea  surface  temperature,  air  temperature  (dry  bulb),  dew- 
point  temperature,  incident  solar  radiatiDn,  and  the  rate  of 
evaporation    (E)    and   precipitatiD  n    (P)  . 

The  turbuien-1-  fluxes  of  latent  heat,  Qe,  and  sensible 
heat,      Qh ,       are    estimated    using    the    bulk    aerodynamic      formu- 
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Qe=   Cd  (.98ES    -    3a)u1Q  (13a) 

Qh=    Cd   (Ts     -    Ta)u10  (13b) 

and    the    net      back   radiation   is    estimated      from   ^.he   empirical 
equation     (Husby    and   Seckel,    1975): 

Qb=    1.14X    10*7(273.  16+TS)4    (.39    -    .5s!/2)  (1    -    .6Q2)  (13c) 

where  Es  is  the  sa-urated  vapor  pressure  of  the  marine  air 
(0.98  corrects  for  salt  defects)  ,  Sa  is  the  vapor  pressure 
of  air  based  on  the  dew-point  temperature,  Ta  air  tempera- 
ture, T3  sea  surface  temperature,  and  Z  is  the  fractional 
cloud   cover.      The   upward    heat    flux,    Qu,    is    then    given    as: 

3U 


Qu=   Qe+    Qh*    2b  (14) 

The    solar    radiation,    Qs,    is  given   by: 

Qg=     (1    -    aocb)  (1    -.66C3)Q0  (15) 

where  QQ  is  the  clear  sky  radiation  given  by  (Seckel  and 
Beauday,    1  973)  : 

Q0=    AQ+    A1  ccs   ^   +    B1  sin  §    +    A2cos    2^    +    B2sin    2^         (16) 

The  constants  a  and  b  are  adopted  from  Tabata  (1964)  and  the 
cubic     cloud   cover      correction   from      Laev=stu    (1960).  The 

variable  is  the  mid-day  elevation  angle  of  the  sun.  The 
coefficients  (A0,A1f  ect . )  were  calculated  by  harmonic  rep- 
resentation of  the  values  predicted  in  the  Smithsonian 
Meteorological   Tables    (list    1953)    with 

^=      (27T/365)   (t-  21)  (17) 

where  t  is  the  Julian  day  of  the  year. 

Not  all  of  the  incoming  short  wave  radiation  penetrates 
the  ocean  mixed  layer.  Approximately  half  (for  open  ocean) 
is  absorbed  within  the  first  leter.  The  amount  absorbed 
varies  from  region  to  region  depending  on  the  amount  of 
absorbing  particles  such  as  ph ytoplankton  and  yellow  sub- 
stance.   More  radiation  will  be  absorbed  in  coastal  regions 
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than  in  the  cpen  ocean  regions.  This  portion  of  absorbed 
short  wave  radiation  is  therefore  considered  to  be  part  of 
the  upward  heat  flux.  The  remainder  of  the  short  wave  radi- 
ation does  penetrate  the  mixed  layer  and  is  attenuated  in  an 
exponential   fashion      depending   on   water   turbidity.  If   the 

fraction  of  solar  radiation  absorbed  in  the  first  meter,  Rf, 
is    given,    then   the   net  heat   at   the   surface    is   given    by: 

Qnet   ■    Qu     *      (Rf)Qs      "    Qs  (18> 

From  these,  the  surface  fluxes  of  buoyancy  (heat  and 
salt)  and  momentum  can  then  be  computed.  The  turbulent  tem- 
perature, salinity,  velocity  and  buoyancy  fluxes  are  given 
by : 

(T'W)  =  -Qnet  <fcp  (19a) 

(S'w«)  =  (P-E)  S0  (19b) 

(u'  w»  )  =  u*  (19c) 

(b'w»)  ~    g[o<(T'w')  -  ^(S'w')]            (19d) 

respectively,  where  u *,  friction  velocity  in  air,  is  calcu- 
lated by: 

rs=j?acdu102  (20) 

u*=     as/pa)1/2  (21) 
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The  subscript  zero  refers  to  the  surface  value  and  ZTS  is  the 
surface   stress       (dynes  cm      ).  In   this    nod  el,        a    positive 

buoyancy  results  when  Qnet<0  and  E>P.  Daring  daytime  peri- 
ods, negative  buoyancy  flux  results  fron  the  domination  of 
the  solar  heating  at  the  surface,  while  during  the  night 
periods,  a  positive  buoyancy  flux  is  produced  by  the  com- 
bined effect  of  net  lcng  wave  radia-ion  and  the  upward  tur- 
bulent  fluxes   of   heat   and    moisture. 

New  ocean  temperature,  salinity  and  wind  driven  current 
profiles    are      predicted   at    one      hour    intervals.  These   are 

then  used  to  predict  the  new  mixed  laysr  dspth.  Ths  steps 
in    the  prediction   computation    are   shown    in    Figure   7. 

D.        FORMULATED    COUPLED   BOUNDARY     LAYER    MODEL 

The  physical  processes  described  in  the  preceding  two 
sections  were  not  altered  during  the  coupling  steps.  Some 
modifications  were  involved  in  transferring  parameters 
between   the   models.  Also   necessary   was    the      uniformity    of 

units.  The  only  parameter  supplied  from  the  ocean  model  for 
use  in  atmospheric  computations  is  the  predicted  sea  surface 
temperature.  However,  this  prediction  is  influenced  by  many 
atmospheric   parameters. 
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Parameters  in  the  atmosphere  initialization  include; 
wind  speed  and  direction,  dew-point  temperature,  latitude, 
Julian  day,  and  time  of  day.  Short  wave  radiative  flux 
(direct  and  diffuse) ,  long  wave  radiative  flux,  sensible  and 
latent  heat  fluxes,  and  the  surface  momentum  flux  were  com- 
puted by  the  atmospheric  section  of  the  coupled  model  and 
input   to    the  ocean   model. 

The  first  problem  associated  with  the  coupling  involved 
the  different  prediction  time  steps  of  the  seperate  models. 
The  atmospheric  model  uses  a  30-minu-e  time  interval, 
whereas,  the  ocean  model  uses  a  one-hour  time  step  interval. 
Assuming  no  significant  sea-surface  temperature  change  would 
rake  place  over  the  short  one-hour  interval,  the  ocean  mod- 
el's prediction  scheme  was  inserted  into  the  atmosphere 
model  as  a  subroutine.  This  prediction  scheme  is  only 
called  every  ether  30-minute  time  step.  The  ccean  model's 
initialization  routines  are  executed  immediately  upen  the 
completion    cf   the  atmosphere   model's    initialization. 

The  second  problem  was  wirh  wind  input.  The  A5L  model 
does  not  require  wind  direction  as  an  input  parame-.er,  only 
ihs  magnitude  of  the  wind  speed.  However,  the  ocean  model 
requires    wind    direction    to   compute   the   horizontal   components 
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cf    the      winds.      These   components      are    used    to      determine  the 

2 

ocean  turbulent  velocity  flux,  u*w  ,  components  used  m  the 
momentum  budget  computations.  A  new  subroutine  was  added  in 
the  initialization  scheme  to  determine  the  horizontal  wind 
components  from  wind  speed  and  direction  input.  These  com- 
ponents   are   determined  by: 

u10x  =  -  (sin  9)u10  (22a) 

uioy  =    "(cos   9)u10  (22b) 

where  u10  is  the  specified  wind  speed  at  10  meters  and  the 
east-west  and  north-south  components  ars  u10x  ar-^  uioy  ' 
respectively.  9  is  the  direction  from  which  the  winds  are 
blowing  relative  to  true  north.  Use  of  aquations  20  and  21 
and   the    following   relationship: 

"Cs=  f^*    =  j>w**w2  (23) 

2 

allows  the  ocean  turbulent  velocity  flux  (u*,  )  tc  be  deter- 
mined (assuming  the  density  of  seawater,  pw  ,  tc  be  approxi- 
mated by  1)  by  first  computing  the  wind  stress  from 
atmospheric  parameters  and  then  equating  the  ocean  turbulent 
velocity    flux    to    the    computed    wind   stress    value.         The    hori- 

2 

zontai  components  of    the    turbulent      velocity    fluxes    i*wx     and 

2 

u*   ,  ar  e  then  easily  computed  by: 
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u*    =Cfi>ii102=    C^(u10x2    +   u10y2)  (2ft) 

u^x2=   C^u10x2  (25a) 

u*wy2=   C$u10y2  (25b) 

The  atmospheric  derived  values  are  computed  using  the  MKS 
unit  system  and  had  to  be  converted  to  tha  CGS  system  prior 
to  any  oceanic  computations.  This  conversion  is  done  at  the 
start  of   each    ocean   prediction   run. 

As  discussed  in  section  B,  the  short  wave  radiative  flux 
is  calculated  using  the  delta-Ed ingt on  method.  The  long  wave 
radiative,  sensible  heat  and  latent  heat  fluxes  are  computed 
using  equations  10a,  10b  and  10c ,  respectively.  A  unit  con- 
version from  Watts  m  to  cal  cm  s  is  required  on  these 
values  prior  to  use  by  ocean  model.  This  conversion  is 
again  done  a-  the  start  of  the  ocean  prediction  run.  A  flow 
diagram  of  the  steps  in  the  coupled  model's  prediction  com- 
putation   is   shown    in    Figure   8. 
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III.       DATA    AND    MODEL    RESULTS 

The  data  sets  used  in  this  study  are  from  the  Coopera- 
tive Experiment  on  West  Coast  Oceanography  and  Meteorology 
from  21  September  to  12  October  1976  (CEW3DM-76)  and  from  5 
to  23  May  1978  (CEWCOM-78).  Both  were  conducted  in  the  Los 
Angeles-San  Diego  basin.  Oceanic  and  atmospheric  data  were 
collected  simultaneously  throughout  th  =  se  time  periods. 
XBT's  were  deployed  from  the  R/V  ACANIA  :>n  a  regular  basis 
(2-3  per  day)  in  conjunction  with  atmosphere  radiosondes. 
Sea  surface  temperature  and  atmospheric  data  such  as  winds, 
temperature,  and  relative  humidity  were  also  obtained  con- 
tinuously   through   these    periods    on   board    the    R/V    ACANIA. 

Three  24-hour  periods  were  selected  for  model  initiali- 
zation and  verification.  The  firs-1:  two  cases  cover  a 
48-hour  period,  5/19/0500  to  5/21/05D0  PST,  during 
CEWCOM-78.  This  period  coinci3es  with  the  observations  dis- 
cussed in  chapter  2,  in  which  coupled  air-sea  interactions 
were  cited  as  possible  explanations  for  ohanges  in  both  the 
atmospheric  and  oceanic  mixed  layer.  Satellite  imagery 
shewed  uniformly  increasing  s-ratus  coverage  (thin  t:o  heavy) 
during     zhis    48-hour      period.         A     clear      sky   situation      was 
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chosen  for  the  third  case.  This  24-hcur  period,  10/04/1900 
to    10/05/1900    PST,    was  during    CEWCOM-76. 

The  general  location  of  the  R/V  ACANIA  during  these 
periods  is  shown  in  Figures  9,  10  and  11.  During  case  1  the 
ship  was  at  anchor  for  the  majority  of  the  time  period.  The 
anchorage  location  is  shown  in  Figure  9.  For  cases  2  and  3 
the  R/V  ACANIA  was  cruising  slowly  (2-3  kiots)  returning  to 
an  initial  point  approximately  every  12  hours.  Winds  were 
primarily  from  the  west  during  the  case  1  period.  The 
winds  shifted  slightly  to  a  west-northwest  direction  midway 
through  the  case  2  period.  Winds  for  case  3  were  primarily 
from  the  northwest  and  were  much  slower  than  in  either  case 
1    or    2. 

Two  predictions  were  analyzed  for  each  case,  one  predic- 
tion was  based  on  the  coupled  air-sea  boundary  layer  model, 
and  the  second  prediction  was  with  the  atmospheric  model 
with   a   fixed    surface    temperature.  A   separate   oceanic   pre- 

diction using  fixed  atmospheric  parameters  (persistence)  was 
not  made.  This  was  the  result  of  an  ir.i_ial  modification 
that  removed  the  wind  stress  and  heat  flux  computation  pack- 
age in  the  ocean  model  in  preparation  foe  coupling  with  the 
atmosphere    model. 
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Figure  9.      Anchorage    Position    of    F/V    ACANTA    in    the   Vicinity 
of    San    Nicolas    Island    from    5/19/78   to    5/20/78. 

Presentation   of  results   will    be    in   two    parts.      The    first 

-co      be  presented      is    the      coupled      model's    performance      with 

respect   to    observed   oceanic   and    atmospheric    variations.       The 

second      will   deal      wi~h      a      comparison    of      predicted      values 

between  the  coupled      model   and   the   atmospheric      model.         All 

rabies   and      figures    are    grouped      together    by    case      number    at 

the    end   of    the   chapter. 
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Figure   10.       Positions   of    P./V    ACANIA   from    5/18/78    to    5/25/'73 
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A.        RESULTS    FROM    COUPLED    MODEL 

1.       CASE    1  (0  500    19   May    £o    0  500    23    Ma£) 

A  listing  of  predicted  atmospheric  and  oceanic  mixed 
layer  quantities  for  this  case  appear  in  Table  I.  Values 
are  for  each  30-minute  time  step  and  since  oceanic  values 
are  only  computed  at  one-hour  intervals,  the  same  value 
appears   at    two   time  steps.  The   parameters   listed   include: 

height  of  inversion  (Zi)  ,  height  cf  lifting  condensation 
level  (Zlcl) ,  wind  speed  (Wind),  atmospheric  mixed  layer 
potential  temperature  (Th) ,  evaporation  duct  height  (Ze) , 
sea  surface  temperature  (SST)  ,  ocean  mixed  layer  depth 
(MLD)  ,    and   the   net   ocean    surface    heat   flux    (TW)  . 

Predicted  values  are  also  sho^n  in  Figures  12 
through  17.  Figure  12  shews  the  time  evolution  of  predic- 
ted sea-surface  temperatures  (solid  line)  and  the  values  of 
the  measured  surface  temperatures  (dashed  line).  Figure  13 
shows  the  time  evolution  of  the  ocean  nixed  layer  depth. 
Figures  14-16  show  ocean  vertical  profiles  for  temperature, 
u  and  v  wind  driven  current  velocities,  and  the  salinity 
change  from  an  initial  value  of  33.5  ppt,  respectively. 
Figure    17    shows    predictions   of   atmospheric    parameters. 
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High  solar  insulation  luring  the  daytime  resulrs  in 
a  net  downward  heat  flux  (ocean  warming)  .  This  can  be  read- 
ily observed  in  the  net  surface  heat  flux  values  (TW)  ,  in 
Table  I,  with  the  negative  sign  appearing  during  the  daytime 
hours.  Maximum  heating  of  tha  ocean  occurred  between  1300 
and  1400  local  time.  This  same  pattern,  warming  in  the  day- 
time hours,  is  also  observed  in  the  sea  surface  temperature 
values,    Table   I   and  Figure   12. 

The  winds  appeared  to  be  fairly  constant,  8ms"1, 
between  0500  and  1100  on  the  nineteenth,  ifind  mixing  coupled 
with  convective  cooling  at  the  surface  caused  the  mixed 
layer  depth  to  deepen  slightly  from  19.5  to  20.9  m  (Figure 
13) .  This  deepening  and  net  heat  loss  at  the  surface  caused 
the  surface  temperature  to  drop.  As  heating  of  the  sea  sur- 
face began  to  increase  due  to  the  increase  of  solar  radia- 
tion,   a  shallow   mixed   layer  of   1 4    meters   formed. 

The  winds  rapidly  decreased  fron  3  to  5  ois"1  at  1200. 
This  coincided  with  a  near  maximum  value  of  radiative  heat- 
ing causing  a  shallow  7  me-er  mixed  laysr  ~o  form  in  less 
than  1  hour.  As  daytime  heating  continued,  a  shallower 
mixed  layer  cf  3.5  meters  formed.  This  process  is  referred 
to    as   the    'afternoon    effect1    in    acoustic    literature   and    is    a 
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quite  common  phenomenon  in  many  areas  of  the  ocean.  In  Fig- 
ure 14,  the  7  meter  mixed  layer  is  barely  discernible  in  the 
1300  vertical  temperature  profile.  The  w  =  ak  gradient  at  the 
base  of  the  mixed  layer  only  required  weak  mixed  layer  sro- 
sion  (increase  in  winds  or  connective  cooling)  to  break  down 
the  shallow  mixed  layer.  During  the  late  afternoon,  the 
winds  began  to  increase,  and  with  the  decreasing  of  solar 
radiation,  the  shallow  thermocline  was  eroded  away.  The  HLD 
continued  deepening  for  the  remainder  of  the  time  period  due 
to    wind- driven   and   convective   mixing. 

Since  sea-surface  temperature  is  the  sole  input  from 
the  ocean  model  to  the  atmosphere  model,  a  comparison  is 
made  of  predicted  and  actual  sea-surface  temperature  values 
(Fiqure  12).  Since  the  R/V  A3ANIA  was  at  anchor  during  the 
majority  cf  this  time  period,  this  is  considered  the  best 
case  to  determine  hov  well  sea  surface  temperature  changes 
due  to  atmospheric  forcing  are  being  predicted.  The  coupled 
model  prediction  yielded  a  0.25 °C  temperature  increase  dur- 
ing the  daytime  heating  period  with  a  maxiaum  temperature  of 
15.4°C  at  1600  local  time.  A  maximum  surface  temperature  of 
15.6°C   was   observed   at   1600   giving   a    0.2°C    error. 
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After  the  daytime  heating  period,  the  surface  tem- 
ps raturss  began  to  drop.  The  minimum  predicted  surface  tem- 
perature at  the  end  cf  the  24-hour  period  was  15.2°C  which 
yields  a  0. 2°C  temperature  drop.  The  0.06°C  temperature 
drop  prior  to  ocean  warming  males  s  the  magnitude  of  the  pre- 
dicted diurnal  warming  and  cooling  to  be  approximately 
egual.  The  actual  minimum  temperature  at  the  end  of  the 
period  was  14.96°C.  Although  there  are  small  errors  in  the 
magnitude  of  the  diurnal  heating  and  cooling  affects,  the 
time    evolution    of  this  phenomenon   is    predicted    quite    well. 

Precipitation  was  not.  observed  during  this  period, 
therefore  the  surface  moisture  flux  and  surface  salinity 
flux  were  due  entirely  to  evaporation.  The  slight  predicted 
increase  in  salinity  from  33.5  to  33.5038  ppt,  shown  in  Fig- 
ure   16,    is   due    to   evaporation. 

An  unstable  condition  existed  in  the  atmosphere  dur- 
ing this  prediction  period.  A  strong  buoyancy  flux  at  the 
surface  increases  entrainment  and  should  oause  the  inversion 
height  to  increase  unless  subsidence  (associated  with  shal- 
lowing of  the  inversion  height)  is  strong  enough  to  prevent 
this  increase.  In  this  case,  subsidence  was  relatively  weak 
(-.0028  ms"  )  and  entrainment  was  able  to  increase  the  inver- 
sion  height. 
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Clcud  occurance/dissipat ion  was  of  interest  because 
cf  its  affect  on  radiation  at  the  surface.  With  the  lifting 
condensation  level,  LCL,  lower  than  the  inversion  throughout 
this  time  period,  stratus  clouds  should  be  present  during 
the  entire  period.  This  feature  was  observed.  During  this 
period  the  LCL  decreased  from  2  43  to  125  meters  while  the 
inversion  height  increased  from  256  to  348  meters.  This 
would  indicate  a  thickening  of  the  clouds  over  this  24-hour 
period  which      was   also      observed    to      be   true.  The   air-sea 

coupling  effect  on  predicting  these  cloud  features  cannot 
be  identified  until  the  results  with  a  constant  sea  surface 
-emperature  forcing  are  examined.  However,  on  the  basis  of 
the  predicted  temperature  change,  a  significant  influence  by 
the    coupling   is   not    expected. 

2.      CASE    2  (0  500    20    May.    to    0  5  00    2_1    May.) 

A  listing  of  predicted  atmospheric  and  oceanic  mixed 
layer  quantities  for  this  period  are  given  in  Table  III  and 
shown  graphically  in  Figures  19  through  24.  As  in  case  1, 
daytime  heating  at  the  surface  is  apparent  by  the  negative 
signs  in  front  of  the  net  surface  heat  flux  values.  The 
maximum  downward  heat  flux  value  was  not  as  high  as  in  case 
1    and   was    probably   due  to   thicker   cloud      conditions    for    this 
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case.  However,        increasing        surface      temperatures      were 

observed    during  J-his    heating   period. 

The  winds  increased  steadily  from  6  to  11  ms" ' 
throughout  the  first  half  of  this  period  and  slowly 
decreased  to  about  9.4  ms"1  at  the  end  of  the  period.  The 
higher  winds  should  be  associated  with  increased  wind  mixing 
cf  the  ocean  mixed  layer.  The  increase!  winds  caused  the 
mixed  layer  depth  to  deepen  at  the  end  of  this  peiod.  This 
also  reduced  the  shallowing  of  the  mixed  layer  during  the 
afternoon  radiative  heating  period.  The  shallowest  mixed 
layer  formed  at  16  meters  versus  3.5  meters  in  the  previous 
case.  In  Figure  21,  the  20  meter  thermooline  is  seen  to  be 
very  weak,  and  with  the  higher  winds,  the  shallower  after- 
noon   mixed    layer   lasted   only   a   couple    of    hours. 

The  effects  of  the  diurnal  heating  and  cooling  peri- 
ods caused  a  0.1°C  temperature  increase  and  a  0.15°C  temper- 
ature decrease,  respectively.  Both  values  are  also  smaller 
than  those  predicted  in  case  1  .  Seasons  for  the  smaller 
predicted  increase  are:  1)  thioker  clouds  reduced  the  amount 
cf  solar  heating  during  the  day,  therby  reducing  the  diurnal 
temperature  increase;  2)  the  higher  winds  increased  the  ver- 
tical  mixing,      causing  cooler   water      below    to    be      mixed    into 
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the  upper  layer  and  also  increased  the  latent  heat  loss  due 
to  evaporation,  thus  reducing  the  total  net  heat  flux  at  the 
surface;  and,  3)  the  larger  air-sea  temperature  difference 
increased  the  sensible  heat  loss,  which  also  reduced  the  net 
heat  flux.  The  decrease  in  the  downward  solar  radiative 
heating  and  increase  in  sensible  and  latent  heat  losses  were 
somewhat  balanced  by  the  decrease  in  the  long  wave  heat  loss 
due  to  the  presence  of  the  cLDuds.  This  decrease  in  the 
long  wave  heat  loss  was  also  the  primary  reason  for  the 
slightly  smaller  value  in  the  nighttime  cooling  period. 

Comparisons  cf  the  predicted  surface  temperature 
values  with  the  actual  values  (Figure  19),  show  a  noticeable 
difference.  This  is  probably  i ue  to  ih£  fact  that  at  0700 
local  time  or.  20  May,  the  R/7  ACANIA  began  cruising  slowly 
to  the  northwest.  The  regional  horizontal  temperature 
structure  must  be  taken  into  account.  The  southern  Califor- 
nia coastal  bathymetry  and  associated  currents  greatly 
influence  the  sea  surface  temperature  structure.  The  Cali- 
fornia Current  brings  in  colder  water  frsm  the  north  as  it 
heads  south  along  the  coast.  Near  Point  Conception,  the 
current  begins  to  turn  to  the  ssuthwest  a*ay  from  'rhe  coast. 
Due  to   coas+al  bathymetry,   Point  Conception   also  provides 
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sheltering  for  the  intrusion  of  warm  water  along  the  coast 
from  the  south  by  the  Davidson  Current.  The  <=ddies  and 
meanders  of  cold  and  warm  water  make  the  sea  surface  temper- 
ature structure  complex  in  the  CEWCOM  experiment  areas.  The 
surface  temperature  variation  observed  by  the  R/V  ACANIA 
appears  to  be  due  to  oceanograph ic  features. 

The  1°C  -"Temperature  increase  in  the  first  two  hours, 
when  the  ship  began  moving  to  the  northwest,  indicates  that 
the  ship's  ^ransit  was  through  one  of  the  small  warm  eddy 
features.  Two  reasons  for  this  are  1)  Radiative  heating  of 
the  ocean  did  not  even  begin  until  after  the  observed  sur- 
face temperature  had  already  begun  to  rapidly  decrease,  and 
2)  the  higher  winds  would  tend  to  mix  the  surface  waters 
with  the  cooler  water  from  below,  thus  causing  a  decrease  in 
temperature  not  an  increase.  The  short  time  duration  of 
this  warm  feature  also  suggests  a  transit  through  an  ano- 
maly. With  the  ship  only  cruising  at  2-3  knots,  the  warm 
anomaly  still  only  lasted  less  than  4  hours. 

The  ship  turned  on  a  northerly  oourse  at  1300  and 
the  rapid  temperature  decrease  stopped.  It  is  reasonably 
safe  to  infer  from  the  sudden  ohange  in  the  surface  tempera- 
ture trend  (rapid  decrease  to  slight  increase)  that  the  ship 
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is  initially  crossing  the  temperature  gradient  contours  and 
then  followed  a  parallel  course  with  these  contours.  This 
would  make  oceanic  conditions  more  compatible  with  those  in 
case    1.  Since  the    0.4°C     temperature   increase      during  the 

daytime  heating  period  compares  favorably  with  the  value 
observed  while  the  ship  was  at  anchor,  the  increase  is 
therefore    presumably    due    to  the    diurnal   heating. 

The  ship  turned  back  to  a  west-northwest  course 
shortly  after  1800  and  the  surface  tempsratur^  immediately 
began  to  decrease  1.5°C  in  4  hours.  Two  hours  after  the  end 
cf  this  period,  0700,  the  ship  came  about  and  paralleled  its 
previous  course  but  in  the  opposite  direction.  An  examina- 
tion of  the  observed  sea  surface  temperatures  shows  a  rapid 
increase  comparable  to  the  decrease  while  heading  in  the 
opposite  direction.  This  further  substantiates  the  observed 
large  temperature  changes  are  not  due  to  atmospheric  forc- 
ing, but  rather  to  the  ship's  movement  through  an  oceano- 
graphic   eddy    or    front. 

Unstable  atmospheric  conditions  existed  during  this 
period  as  well.  The  1.5°C  drop  in  the  air  temperature  and  a 
much  smaller  drop  in  the  ocean  temperature  (0.2°C)  coupled 
with    higher   wind   speeds   resulted    in    a    large    surface   buoyancy 
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flux.  This  resulted  in  a  much  higher  predicted  inversion 
height  change  due  to  en t raiment .  The  slightly  smaller  sub- 
sidence rate  was  also  a  factor  in  the  greater  predicted 
increase    of   the    inversion    height. 

The  lifting  condensation  level  fluctuates  up  and 
down  unlike  the  steady  decrease  in  height  in  c*se  1.  How- 
ever, the  height  of  the  LCL  never  increased  above  the  ini- 
tial height  and  so  with  the  large  increase  in  height  of  the 
inversion,  thickening  of  the  stratus  cloui  deck  is  predicted 
to    continue  during  this   period   as   wall. 

3.       CASE    3  (JiQO    04   OCT    to    1900    05    OCT) 

A  listing  of  predicted  atmospheric:  and  oceanic  mixed 
lay^r  guantities  are  given  in  Table  V  and  shown  graphically 
in  Figures  26  through  31.  Although  the  time  period  for  this 
case  is  in  October,  there  are  still  many  features  occuring 
that  are  similar  to  the  first  two  cases  ii  May.  The  primary 
differences  in  this  case,  other  than  being  a  cloud-free 
case,  are  1)  the  air  and  surfaos  temperature  are  much  warmer 
and    2)    a    much    higher    value   of    subsidence    was    observel. 

Daytime  heating  is  apparent  in  predictions  for  the 
clear  sky  case,  as  expected.  The  magnitude  of  J-he  net  ocean 
surface    heat    flux    is,      however,         smaller    than    the    value    for 
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the  thin  stratus  deck  condition  in  case  1.  This  is  presum- 
ably due  to  the  difference  in  seasonal  heating  rate  and  will 
not    be  considered   part  of   the    local   coupling   effects. 

The  winds  in  case  3  are  much  weaker  than  in  the 
first  two  cases.  With  the  Light  winds  and  high  downward 
surface  heat  flux  during  the  day,  a  much  shallower  mixed 
layer  might  be  expected.  The  ocean  mixed  layer  depth  (MLD) 
values  in  Table  V  (also  shown  graphically  in  Figure  25)  show 
that  an  ocean  mixed  layer  of  1  meter  suddenly  forms  in  less 
fhan    an      hour   after      heating  began      at   0930.  This    shallow 

thermocline  is  much  stronger  than  in  either  case  1  or  2,  and 
persists    throughout      most    of    -tie      day    (Figure    28).  As    the 

winds  slowly  increased  late  in  the  afternoon,  the  shallow 
mixed  layer  deepens.  Without  winds  to  mix  this  shallow 
region,  surface  heating  is  much  greater  in  this  case.  A 
0.6°C  surface   temperature    increase   was   predicted    to    occur. 

As  in  case  2,  the  ship's  movement  played  a  major 
role  in  -che  large  variations  in  observed  surface  temperature 
values.  This  prevented  the  verification  of  the  predicted 
surface  temperature  increase  during  the  daytime  heading 
period.  Figure  11  shows  the  complex  horizontal  temperature 
structure    in      the  CSWCOM-76      experiment   area      for  the      first 
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week  of  October.  From  this  figure,  it  is  hard  to  explain 
the  large  2.8°C  temperature  drop  during  the  first  half  of 
the  period  as  the  ship  appeared  to  be  heading  into  warmer 
water.  The  rapid  temperature  increase,  starting  at  0900,  is 
obviously  due  to  the  ship  transiting  into  the  warm  eddy  fea- 
ture  north   of   Santa   Catalina   Island    {Figure    11). 

Unstable  atmospheric  conditions  are  also  present 
during  this  time  period.  The  high  subsidence  caused  the 
inversion  to  be  forced  downward  to  within  72  meters  of  the 
surface.  The  lifting  condensation  level  decreases  even  more 
rapidly  and  finally  becomes  lower  than  the  inversion  at  the 
last  hour  of  the  forecast  period.  This  implies  the  forma- 
tion of  a  low  stratus  cloud  deck.  From  Figure  33,  fog  and 
stratus  clouds  were  observed  to  form  in  the  CEWCOM  experi- 
ment   area. 

E.        COMPARISON    OF    COUPLED    AND    AT  MOSPHERIC    MODEL 

The  second  and  final  phase  of  this  study  is  to  compare 
the  previous  coupled  air-sea  boundary  layer  model  results 
with  the  atmospheric  model  results  using  a  constant  sea  sur- 
face temperature.  This  comparison  is  to  determine  if  any 
differences  existed  between  the  atmospheric  predictions,  and 
if    so,        determine   if  these      differences    are    due      to    coupled 
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interactions.  A  listing  of  the  atmospheric  model  predic- 
tions of  atmosphere  mixed  layer  parameters  for  cases  1,  2 
and  3  are  found  in  Tables  II,  IV  and  VI,  respectively. 
Graphical  presentations  of  these  results  are  shown  in  Fig- 
ures 18,  25  and  32.  All  tables  and  figures  are  grouped 
together    by   case   number    at   the    and   of  this    chapter. 

The  initial  comparison  reveals  no  significant  difference 
in  the  two  predicted  sets  of  mixed  layer  parameters,  indi- 
cating a  negligible  effect  of  coupled  interactions  on  atmos- 
pheric predictions  for  xhis  short  time  period.  The 
predicted  inversion  heights  are  essentially  identical  in 
both  sets  of  predictions  for  all  three  oases  (maximum  dif- 
ference of  2  meters) .  This  is  also  found  to  be  true  for 
humidity  and  even  the  evaporation  duct  height  predictions. 
This  observation  is  somewhat  surprising  since  there  are  sig- 
nificant interactions  between  the  surfaoe  fluxes  (evapora- 
tion) and  these  parameters.  The  close  proximity  of  these 
parameters  to  the  ocean  surface  would  also  imply  a  more  sig- 
nificant impac"-.  even  with  small  changes  in  sea-surface  tem- 
peratures. 

The  only  parameter  than  seems  to  be  olearly  affected  by 
the    coupling      is   the    lifting   condensation      level.         However, 
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even  for  the  cloudy  conditions  of  cases  1  and  2,  this  shows 
only  a  slight  deviation.  A  maximum  difference  of  only  5 
meters  is  observed  between  the  predictions  for  these  cases. 
In  the  clear  sky  case  under  warmer  and  more  humid  condi- 
tions, there  is  a  30  meter  difference  observed  in  the  two 
LCL    predictions.  However  this    occurred    as      the   atmosphere 

approach    near      neutral  conditions.  This    may      indicate   the 

model's  inability  to  make  reliable  forecasts  under  stable 
atmospheric  conditions. 

Since  a  comparison  of  predicted  ocean  parameters  between 
a  coupled  and  independent  ocean  model  was  net  made,  no  con- 
clusions on  this  topic  can  be  made  at  this  time.  Although 
tentative  results  indicate  no  definite  advantages  using  a 
coupled  model  for  the  atmospheric  predictions,  the  reverse 
may    very    well   be   true    for    oceanic   predictions. 
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IV.       COLLUSIONS/RECOMMENDATIONS 

On  the  basis  of  a  limited  study  of  several  24-hour  pre- 
dictions, surface  temperature  variations  are  not  large 
enouqh  to  require  ocean  mixed  layer  processes  in  short  term 
atmospheric   forecasts.  Carefully   selected      data   sets      and 

considerable  interpretive  efforts  are  essential  if  general 
conclusions  are  to  be  drawn  on  the  relative  roles  of  dynamic 
processes      in    bo^h      boundary   layers.  Further   studies      are 

required  to  determine  the  model's  ability  to  better  handle 
forecasts  at  longer  time  scales.  Further  studies  should 
also  include  sensitivity  analyses  to  define  situations  for 
which  air-sea  interaction  effects  are  greatest.  Studies  of 
oceanic  sensitivity  to  atmospneric  forcing  are  required  to 
determine    the    impact    it    will   have   en    acoustical    forecasts. 

The  CEWCOM-76  and  78  data,  although  able  to  supply 
simultaneous  and  oceanic  lata,  suffered  in  two  major 
respects  First,  the  complexity  of  the  horizontal  temperature 
structure  with  meanders  cf  wan  and  cold  eddy  features  had 
definite    impacts   on   observed    data.  Second,      data    was    col- 

lected while   the    R/V    ACANIA  was    transiting    through   a   complex 
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ocean  temperature  structure.  Changes  in  observed  conditions 
under  these  circumstances  can  and  did  mask  the  interaction 
processes   of   interest. 

Although  the  primary  results  of  this  study  indicate 
coupled  air-sea  models  are  probably  not  necessary  for  short 
term  atmospheric  forecasts,  a  coupled  modal  still  has  advan- 
tages, especially  if  long  term  forecasts  are  to  become  pos- 
sible. Such  a  model  would  not  be  the  single  one-dimensional 
one  considered  here.  However,  this  model  offers  many  poten- 
tial benefits  as  a  multi-purpose  model  with  single  station 
assessment  and  prediction  capabilities  in  the  absence  of 
shore-based  support.  A  sound  vsloci-y  profile  (SVP)  predic- 
tion capability  could  be  easily  incorporated  into  this 
model.  In  conjunction  with  present  prediction  capabilities, 
a  multi-purpose  model  would  be  available  for  a  wide  range  of 
needs.  Ocean  mixed  layer  depch  and  SVP  predictions  can  be 
used  in  predicting  ASW  sensor  oapabiliti? s .  Environmental 
effects  on  the  performance  of  weapons  and  sensors  that 
depend  on  electromagnetic  and  electro-optical  wave  propaga- 
tion   can    be   evalua-ed   using  the    same    modal. 
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Figure    12.      Time    Evolution   of    the    Sea-Surface    Temperature- 
Case    1.       Predicted   values    are    shown    by   the    solid 
line.      Observed  values   are    shown    by   the   dashed 
line. 
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Figure  13.   Time  Evolution  of  Ocean  Mixed  Layer  Deoth  (MLD)  - 
Case  1 . 
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Figure   14.      Selectee1.    Ocean   Vertical   Temperature    Profiles- 
Case    1 . 
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Figure    15.       Selected    Ocean   Vertical   Profiles    of    the 

Horizontal  Wind   Driven   Currents-    Case    1.      Fas- 
west    component   u    (solid    line)       North-south 
component   v    (dashed)  . 
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Figure  17.   Graphical  DisDlay  of  Lh=  Coupled  Model  Predicted 
Atmospheric  Mixed  Layer  Quantities-  Case  1. 
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Figure   18.      Graphical  .Display   of    the    MA3L    Model    Predicted 
Atmospheric    Mixed    Layer   Quantities-    Case    1. 
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Figure   19.      Time   Evolution  of   the   Ssa-Surface   Temperature- 

Case   2.       Predicts*!   values    are    shown    by   the    solid 
line.      Observed  values   are   shown    by   the   dashed 
line. 
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Fiaure  20.   Time  Evolution  of  Dcean  Mixed  Layer  DeDth  (MLD) 
Case  2. 
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Figure   2  1.       Selected    Ocean   Vertical  Temperature    Profiles- 
Case    2. 
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Figure   22. 


Selected    Ocean   Vertical   Profiles    of   the 
Horizontal  Wind   Driven   Currants-    Case    2.       Eas' 
west    component   u    (solid    line)       North-south 
component   v    (dashed)  . 
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Figure  24.      Graphical  .Display   of    the    Coupled    Model   *r-d-c^d 
A-mospheric    Mixed   Layer   Quantities-   Case    2.       "~ 
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Figure  25.   Graphical  Display  of    ■uh=  MABL  Model  Predicted 
Atmospheric  Mixed  Layer  Quantities-  Case  2. 
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Figure   26.      Time    Evolution  of    the   Sea-Surface   Temperature- 
Case   3.       Predicted   values    are    shown    by   "he    solid 
line.      Observed  values   are   shown    by   tae    dashed 
line. 
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Figure   27.       Time   Evolution  of    3cean    Mixed    Laver    Depth     (HLD) - 
Case    3. 
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Figure   28.      Selected    Ocean   Vertical  Temperature    Profiles- 
Case   3. 
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Figure   29. 


Selected    Ocean   Vertical   Profiles    of    the 
Horizontal  Wind   Driven    Currents-    Case    3.       Fas' 
west    component   u    (solid    line)       North-south 
component   v    (dashed)  . 
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Fiqure   30.       Selected    Change   in   Dceaa    Vertical    Salinity 
Profiles-   Case   3. 
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Figure  3  1.   Graphical  Display  of    +  he  Couolsd  Model  Predicted 
Atmospheric  Mixed  Layer  Quantities-  Cass  3. 
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Figure  32.   Graphical  Display  3f  the  MABL  Model  Predicted 
Atmospheric  Mixea  Layer  Quantities-  Case  3. 
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